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We study the melting of charge order in the half doped manganites using a model that incorporates 
double exchange, antiferromagnetic superexchange, and Jahn- Teller coupling between electrons and 
phonons. We primarily use a real space Monte Carlo technique to study the phase diagram in terms 
of applied field {h) and temperature (T), exploring the melting of charge order with increasing h and 
its recovery on decreasing h. We observe hysteresis in this response, discover that the 'field melted' 
high conductance state can be spatially inhomogeneous even in the absence of extrinsic disorder and 
find that the hysteretic response plays out in the background of a curious equilibrium physics of 
the melting transition. Our extensive numerical work, exploring h, T, and the electronic parameter 
space, is backed up by discussion of simpler limiting cases, and a Landau framework for explaining 
the hysteresis and field induced transitions. We compare our results in detail to experimental data 
on the half doped manganites. This paper focuses on our results in the 'clean' systems, a companion 
paper studies the effect of cation disorder on the melting phenomena. 



I. INTRODUCTION 

The manganites, LnMnO^, where Ln belongs to the 
lanthanide family, display a rich variety of phase^^ 
as a function of hole doping and bandwidth. Promi- 
nent among these are the ferromagnetic metal (FM- 
M) and the antiferromagnetic charge ordered insulat- 
ing (AF-CO-I) states. The ferromagnetic metal usually 
shows up in manganites with large bandwidth, i.e, large 
mean cation radius, rA, for hole doping x ~ 0.2 — 0.5, 
while the AF-CO-I state is observed in low bandwidth 
materials at commensurate doping, x ~ 0.5. This 
'half doped' state is especially interesting since it al- 
lows systematic study of phase competition and the 
role of disorder. The x = 0.5 state ha s been exten- 
sively probed expcrimcntall^l^llIIIIEIIIIIini and also anal- 
ysed theoretically" ^offl 

It is known that low bandwidth (BW) materials have 
CE magnetic order with a related checkerboard charge 
order and concomitant orbital order are insulating (CE- 
CO-I), while the large BW materials have a FM-M 
ground state. At intermediate BW some materials have 
'A type' magnetic order. Application of a magnetic field 
can melt the charge order and convert the CE-CO-I to a 
ferromagnetic metal. The melting transition appears to 
be abrupt an d is accompanied by hysteresis in response 
to field swee^^. 

The general understanding of field induced melting of 
charge order is that an applied field destabilises the CE 
magnetic order in favour of ferromagnetism, and the in- 
timate coupling of charge order and CE magnetic order 
results in the 'melting' of charge order. Some aspects 
of the melting problem are well studied, (i) The ther- 
modynamic melting field he (which is bracketed by the 
actual switching fields hf discussed later) is small. The 
associated energy scale g^-shc ^ kBTco, where Tco is 
the zero field melting temperature. The smallness of he 
is attributed to the small energy difference between the 
CE-CO-I and FM-M states, (ii) The field induced tran- 



sition is believed to be first order in nature and is accom- 
panied by hysteresis arising out of competing metastable 
minima, particularly for materials close to the CE-CO-I - 
FM-M phase boundary, (iii) The field induced FM-M is 
believed to be homogeneous, and is so assumed in theo- 
retical studies. These general observations set important 
benchmarks for any detailed theory, but some key issues 
remained unresolved. 

(i) The nature of the finite field state: While the CE- 
CO-I state melts on applying a magnetic field, it was 
not clear that the g round state is a homogeneous metal. 
Recent experimentJ^l^tSIZI jn fact demonstrate that the 
melted state in actually inhomogeneous. We explicitly 
show that at half doping, for bandwidth materials with 
weak CE-CO-I state, phase separation at intermediate 
fields is inevitable and the 'melted' state is at best a 
percolative metal. It is only at large fields, does one 
recover the homogeneous FM-M final state. Thus, con- 
trary to the belief that the field melting is a discontinu- 
ous transition, we contend that the transition is actually 
continuous, with equilibrium phase separated states in- 
terpolating between the initial(/i = 0) CE-CO-I and the 
final homogeneous FM-M states. We suggest that in real 
materials with low BW and very weak disordei^^^, the 
continuous nature of the transition could be probed. 

(ii) Switching fields and limits of metastability: While 
the thermodynamic critical field can be estimated from 
energy balance, the upper and lower critical fields that 
are actually measured define limits of metastability, and 
have remained inaccessible within the calculations per- 
formed so far. 

(iii) The effect of disorder on field melting: One ex- 
pects he to increase with reducing bandwidth, since CO 
is better stabilised. This indeed happens for lanthanides 
(Ln) of the forrrPLni/2Cai/2Mn03. For members of the 
Ln]^/2Sri/2Mn03 familyl^, with very similar bandwidths, 
he increases initially with decreasing BW but takes a 
downturn beyond a critical BW and then drops to zero, 
in sharp contrast to the 'divergence' of he seen in the Ca 
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FIG. 1: Colour online: The h — T phase diagram of various REi/2AEi/2Mn03 compounds^. The materials involve a systematic 
decrease in ta from Pro.sSro.sMnOa (PSMO) to Smo.sCao.sMnOs (SCMO). The critical CO mehing temperatures increase 
with decreasing rA- The associated hysteresis with the first order nature of transition, opens a window of metastabihty at 
low T, which tapers with increase in temperature and vanishes at Tco- The Ca family has low disorder and shows re-entrant 
behaviour in /ij, which vanishes at very low bandwidths (Smo.sCao.sMnOs). Further, in the Ca family, the decrease in rA, also 
makes the CO state more robust, with SCMO having the largest melting fields. The Sr family has larger disorder (see text) 
with Pro.sSro.sMnOa having larger rA than Ndo.5Sro.5Mn03(NSMO). While NSMO shows marked hysteresis, PSMO is rather 
benign. 



family. We provide an explanatioiJ^ for this in terms of 
the reduced stiffness of the CO state in the presence of 
disorder, and consequently its melting at a lower field. 

In the present work, we have used a real space Monte 
Carlo techniqu^^ to study the melting of charge order. 
In our earlier worlP^, we have addressed some of the 
above mentioned issues, however, this is a much detailed 
exposition of the phenomenon of field melting of CO in 
the half doped manganites. We have mapped out the 
h — T phase diagrams capturing both hysteresis and re- 
entrant features as seen in experiments. In addition to 
thermodynamic indicators we have characterized the sys- 
tem at low temperature through direct spatial snapshots 
and by measuring the volume fraction of the CO-I and the 
FM-M. Using these tools we have explored the combined 
effects of bandwidth variation and disorder on the process 
of field induced melting. We have addressed all the unre- 
solved issues, (i)-(iii), listed earlier. Further, by putting 
our results within a Landau-like energy landscape, we 
have provided a broad framework for rationalizing the 
observed material systematics. 

The paper is organised as follows. In section II we sum- 
marize the key experimental results, and follow it in sec- 
tion III with a discussion of earlier theoretical work on 
field melting. In section IV we define our the model and 
describe the method for solving it. Section V discusses 



the zero field reference state. Section VI is the heart of 
the paper and discusses the results at finite field. The 
later sections, VII and VIII, looks at these results from 
alternate simpler calculations, puts forward a Landau en- 
ergy landscape, and discusses the robustness of our nu- 
merical results. Section IX concludes the paper. 



II. EXPERIMENTAL RESULTS 



Ca Sr Ba 
1.34 1.44 1.61 

La Pr Nd Sm Eu Gd Tb Ho Y 
1.36 1.29 1.27 1.24 1.23 1.21 1.20 1.18 1.18 



TABLE I: Ionic radii (in Angstroms) for various k¥?+ and 
RE''"'' ions in the perovskite manganites^. 

Experimental investigations on the field melting of the 
CO state has probed the hysteretic melting, the band- 
width dependence of the critical melting fields and sys- 
tematic effects of disorder, primarily through transport 
measurements and direct spatial imaging of the lattice 
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undergoing the first order melting transition. We briefly 
review the important results here. The key experimen- 
tal magnetic field-temperature {h — T) phase diagrams 
are shown in Fig.l. The material systematics were ob- 
tained by systematic decrease in ta from Pro.sSro.sMnOa 
(PSMO) to Smo.sCao.sMnOg (SCMO). The phase dia- 
grams were obtained by sweeping up and down in mag- 
netic fields, at fixed temperatures, on samples that have 
been initially zero field cooled. Fig.l shows that the melt- 
ing field in the upward sweep {h^) differ the most from 
the 'recovery' field {h~) in the downward sweep when 
T — !■ 0. This window narrows and vanishes as T — ^ Tco- 
The evolution of h~ has a 're-entrant' feature that van- 
ishes for small BW materials and is pronounced for the 
intermediate and low BW. Further, reducing BW (or re- 
ducing Ta) progressively increases the stability of the CO 
state, with SCMO having the largest and Tco- 

For the experimental results shown in Fig.l, two 
lanthanide families have been extensively investigated, 
the Ca serie^, Lno.sCao.sMnOa , and the Sr serie^, 
Lno.5Sro.5Mn03. The bandwidth is varied by making 
materials with Ln atoms with progressively smaller radii 
(see Table 1). Changing the Ln atom also changes the 
cation size mismatch a a, varying the structural disorder. 
The disorder arises because random substitution of Ln by 
Ca/Sr leads to random changes in the Mn-O-Mn bond 
angles, 6 say, modulating the local hopping, cx cos^{9). 
The other effect of random substitution is that if the net 
charge on the dopants is different than that on the parent 
Ln, then charge disorder creates scattering potentials. 

Traditionallj'Sl the extent of structural disorder is 
quantified by the variance a a of the average ionic radii 
of the A site in Ai-^jA'^^BOs and h as been measured for 
the Ca, Sr, and Ba familie^^^'i^ at various hole dop- 
ings close to half filling, x = 0.45. While in the Ca 
family a a is ^ 10^"^ A^, for the Sr family, greater size 
mismatch (Table. 1) leads to a a ^ 10~^A^. In our work, 
we will consider the Ca family to represent the 'clean' 
manganite, while the Sr family is typical of a disordered 
manganite. 

In addition to BW variation, explored in detail, as for 
example in Fig.l, there are studies that focus on the role 
of disorder in field induced melting. The impact of dis- 
order on the zero field a; = 0.5 state has been beautifully 
demonstrated^. Experiments at half doping where CO 
melts on applying a magnetic field has been carried out 
with a particular focus on the spatial nature of the melt- 
ing proces^^- and the role of disordei^^. 

While in the Ca family the melting fields increase with 
decreasing r^, in the Sr family they are suppressed on 
decreasing r^. For the Br family, with the largest disor- 
der among the three families, the global CO correlations 
are destroyed even in the zero field ground state. Fur- 
ther, for the former two families, the systematic change 
in the upper and lower critical fields with rA, provides 
information on the evolution of the various competing 
phases as the bandwidth changes. The phase diagrams 
in Fig.l, however do not provide details of the spatial na- 



ture of the field melted state. Also the detailed thermal 
evolution, in particular, the high temperature phases are 
not explored. 

Some recent experiments have begun exploring the na- 
ture of the field melted state. The high field state in half 
doped NSMO has been probed for spatial inhomogeneity 
of the magnetic order in the ground stat^. This study 
reveals that the low T, the high field state is rather inho- 
mogeneous with 'poor FM' domains coexisting with per- 
fect FM regions. Since resistivity data implies this state 
is metallic, the system appears to be a percolative metal. 
Similar results were reported in La5_„Pr,,Ca3 MnO-j 

8 y " 8 

where again CO-I regions are shown to coexist with FM- 
M regions. Further experiment^ track the evolution of 
the charge ordered state with magnetic field and temper- 
ature in LCMO, PCMO, PSMO and NSMO at a; = 0.5. 
The main body of results indicate coexistence, at low 
T, of competing (AF-I/FM-M) phases. Further, it is re- 
ported that the fraction of the two phases can be tuned 
by following different protocols in /i — T variation. The 
tunability is ascribed to the disorder induced arrest of 
kinetics of first order phase transition, and is a non- 
equilibrium effect. 

While we address the results on the 'clean' systems 
such as Ca family in this work, we contrast these results 
with that in the disordered systems in the companion 
papeilSl. 



III. EARLIER WORK 

These intriguing experimental results have led to 
attempts at a theory of the field melting in CO 
manganite!|i™Mni2ll 

(i) The earliest attemplfi^involved the mean field study 
of a one band model with on-site and nearest neighbour 
Hubbard interaction in addition to double exchange. On 
application of a magnetic field the zero field AF-CO-I 
state was shown to melt to a FM-M through a first order 
transition. The result shows the existence of a magnetic 
order dependent CO state, which is destabilised when the 
AF order is disrupted by a field. 

(ii) A more elaborate and realistic modeP^l incorpo- 
rating Jahn- Teller interaction has also been explored. A 
variety of magnetically ordered states (FM, G-AF, CE- 
AF) along with a metallic or a CO state, were examined 
in the presence of a field. The destabilisation of the AF- 
CO state with increasing h was tracked for varying BW, 
and an increase in the field needed for the transition with 
decreasing rA was captured. 

(iii) Finally, a two orbital model was studied, with a 
large family of variational states in a recent work^^. This 
established that the smallness of the (thermodynamic) 
melting fields is due to the closeness in energy of the CE- 
CO-I and the FM-M phases. It also suggested that the 
field melted state could be inhomogeneous, since, since 
for range of electron-phonon couplings, at finite h, a CE- 
CO state with defects seems preferable to the homoge- 
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neous CE-CO-I or FM-M. We believe that the result hints 
at field induced phase separation but the authors did not 
pursue the issue further. 

The attempts above have added valuable insight to the 
field melting problem but have left three issues unex- 
plored: (i) the spatial character of the high field state, 

(ii) the impact of disorder on the melting process, and 

(iii) the dependence of the switching field hf on band- 
width and disorder. We use a real space Monte Carlo 
method on large lattices to settle these questions. 

IV. MODEL AND METHOD 
A. Model 

For our calculations, we consider a two band model 
for Cg electrons, Hunds coupled to t2g derived core spins, 
in a two dimensional square (2D) lattice. The electrons 
are also coupled to Jahn- Teller phonons, while the core 
spins have an AF superexchange coupling between them. 
These ingredients are all necessary to obtain a CE-CO-I 
phase. 

{ij)a i {ij) 

~A^Q,.x,, + |^Q2-^7V-;i^5,, (1) 

i i i 

Here, c and are annihilation and creation operators 
for Cg electrons and a, (3 are the two Mn-e^ orbitals 
dx'^-y2 and d^^^-r^, labelled (a) and (&) in what follows. 

hopping amplitudes between nearest-neighbor 
sites with the symmetry dictated form: t^^ — t^^ = t, 

where x and y are spatial directions The eg electron spin 
is erf = X]"<T' c\ac^'^„,Ciaa', whcrc the F's are Pauh ma- 
trices. It is coupled to the t2g spin S,; via the Hund's 
coupling Jh, and we assume Jn/t ^1- A is the cou- 
pling between the JT distortion Q; — {Qix, Qiz) and the 

orbital pseudospin rf = J2a^ '^laa^ap^iP'y^ ™d K is the 
lattice stiffness, and h the magnetic field. We assume 
it to be in the z direction and coupled only to S^. We 
set i = 1, if = 1, and treat the and Si as classi- 
cal variables. The chemical potential /i is adjusted so 
that the electron density remains n — \/2 which is also 
x=l-n = 1/2. 

B. Parameter space 

Even in the model without disorder, the parameter 
space with Jh, J and A is rather large. So one looks to 
restrict the same. In this quest, the large Hund's rule 
coupling, as seen in experiments is assumed to be large 
enough to completely enslave the electron spin to the 



core spin at any site. In the context of the model, this 
amounts to working in the Ja/i oo limit, and thus the 
relevant parameters that remain are the electron phonon 
coupling (A) and the AF superexchange scale (J). In the 
present paper we do not consider disorder so the only 
electronic parameters are A and J. We choose A and J 
such that the system is a CE-CO-I at low T and h = 0, 
but close to a FM-M phase. Changing A is equivalent 
to change in (inverse) BW. This variation is needed to 
model the materials systematics with changing rA- In 
reality J also changes with changing BW but we will 
assume these two to be independent parameters and will 
explore only a couple of J values, J = 0.10, 0.12, to 
narrow the parameter space. The justification of these 
values will be provided in the text. Disorder introduces 
another 'dimension' to the parameter space, those results 
are reported separately. 



C. Method 

The present work is in the spirit of many other 
calculation^^H^ where classical Monte Carlo technique has 
been employed to explore the nonequilibrium effects like 
hysteresis etc., around a phase transition. As mentioned 
earlier, this study aims at describing the 'real' space na- 
ture of field melting and to that end we a real space ex- 
act diagonalization (ED) based Monte Carlo (MC) that 
tracks the spin-charge-orbital degrees of freedom of the 
system across the melting transition. This technique de- 
fines a region (cluster) around the site at which one at- 
tempts a MC update and accepts and rejects the update 
based on the energy change only in the cluster. This 
'travelling cluster approximation' (TCA^^ gives a huge 
gain compared to standard ED-MC the cost for which 
scales as N^, N being the size of the system. Within 
TCA the cost for the same system is ^ NN^, where Nc 
is the fixed cluster size, and is linear in iV as opposed to 
N'^. Using this technique we have accessed sizes up to 
40^ as opposed to the limit of 8^ within ED-MC. In 
the present work we employ 8^ moving cluster on system 
sizes between 16^ (for most data presented here) to 40^ 
(for detailed study of spatial evolution). For calculation 
of observables, we diagonalizes the full system, after the 
system has equlibriated. In studies such as this, where 
real space phase separation may play vital role in de- 
termining transport responses, it is pertinent to be able 
to access system sizes that are large enough to capture 
coexisting phases on the lattice. Such a study is not pos- 
sible reliably within conventional ED-MC due to severe 
size limitations. Most results are obtained by sweeping 
the system in magnetic fields at low temperature. The 
issues of comparison between experimental and Monte 
Carlo sweep rates, etc., has been discussed later in the 
paper. 
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FIG. 2: Colour online: (a) The A — J phase diagram at T = 0.01. The various phases are indicated in colors separated by 
solid lines. The dashed/dotted lines, demarcate various parameter regimes in the CE-CO-I phase terms of their response to 
the applied fields. Clearly, the A, J window of interest lies where the FM-M and the CE-CO-I phases are close by, this happens 
for A ~ 1.6 and J ~ 0.1. For large A values, above the big-dashed line, the CO state doesnot melt on applying a magnetic field. 
Below this line, the CE-CO-I state melts in response to an applied magnetic field. However, this region can be classified into 
one that recovers the CO state in a low T field cycle (entire CE-CO-I region between the big-dashed line and finely-dashed 
line) and another that does not recover CO (the CE-CO-I region below the finely-dashed line). Below the dotted line the field 
melting yields a homogeneous FM-M and above it, the melting forms an inhomogeneous state, (b) The clean A — T phase 
diagram, at J = 0.1, and (c) for J = 0.12. Note, in panels (b) and (c) A increases from right to left. Since, these are thermal 
evolution with A at the fixed J, the AF(A-2D) region separating the FM-M and CE-CO-I phases at low T, is smaller for the 
J = 0.1 cross section as compared to that at J = 0.12. 



D. Physical quantities 

In order to study the evolution with applied magnetic 
field, we track various physical quantities in real space 
and momentum space which provide information about 
the correlations of (and between) the various d.o.f. of 
the system. These are further backed up with real space 
snapshots of the lattice across the melting transition. 
Here we briefly discuss these indicators. 

We compute the 'one point' distribution of lattice dis- 
tortions, P{Q) = J2i^iQ - Qi)^ where Qi = \Qi\, spa- 
tial Q-Q correlations, Dgiq) = E^J(QiQJ)e'^■^""'"""'^ 
and spin-spin correlations, S{q) = ^j^.(Si.Sj)e**^'^'"'~'"j'. 
Angular brackets represent a thermal average. We also 
compute the volume fraction of the charge ordered re- 
gion in the lattice from direct spatial snapshots of the 
charge distribution. To measure the volume fraction, we 
tag a site with a particular color if the site has n > 0.5 
and is surrounded by the four nearest neighbor sites with 
n < 0.5 and vice verse (i.e. a site with local anti-ferro- 
charge correlation is marked with a particular color). 
Similarly, if the difference between the charge density at 
a site with its nearest neighbours is less than a threshold, 
that site is tagged by a different color, i. e. , the charge uni- 
form regions are marked by this color. For intermediate 
cases, we use an interpolative colour scheme. A measure 
of the volume fraction is necessary for studying inhomo- 
geneous melting because the momentum space structure 
factors are not a good measure of the local CO in the 
system. Further, the spatial snapshots also directly pro- 



vide visual information on the melting process. While 
the indicators above measure the correlations and spatial 
evolution, the metallic or insulating character is tracked 
via (low frequency) conductivitjPSl^ Udc, and the density 
of states (DOS), N{uj) — {^J2n^(^~^n)), where e„ are 
the electronic eigenvalues in some MC background and 
the angular brackets indicate thermal average. We track 
all the above quantities as a function of temperature and 
applied magnetic fields for studying the CO melting phe- 
nomenon. 



V. RESULTS AT ZERO FIELD 

Prior to studying the effects of magnetic fields on the 
CE-CO-1 state, it is useful to look at the zero field sit- 
uation. This will set up a reference state, give informa- 
tion about the phases that compete with the CE-CO-1 
phase and help us to justify an appropriate A, J parame- 
ter regime for calculations presented later. Accordingly, 
in this section we discuss in detail, the A — J phase dia- 
gram at low T and h — 0. Further, this phase diagram 
with all the zero field competing phases also serves as a 
plot on which one can not only classify parameter regimes 
with varied responses to a typical magnetic field sweep, 
but also rationalize them to a large extent purely in terms 
of proximity to zero field phase boundaries. 
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A. A — J phase diagram at T = 

Fig. 2. (a) shows the A — J phase diagrairP^ where the 
various phases obtained are indicated. These first or- 
der phase boundaries (indicated by sohd hues) are de- 
termined by cooHng the system in zero field. The 
dashed/dotted hues, obtained by low temperature field 
sweeps, indicate the boundaries of the regions showing 
different responses to field sweeps. While the detailed 
indicators allowing this classification (of regions) will be 
the subject of the next section, here we will just present 
the classification and rationalize it, to the extent possi- 
ble, purely in terms of competing phases and proximity 
to phase boundaries. 

Let us begin by describing the various competing 
phases. The CE-CO-I phase and the FM-M phases 
are separated by a thin region (A-2D) with 'line-like', 
q — {0, tt} spin correlations (Fig. 2. (a)). If we are at a 
parameter regime where the FM-M and CE-CO-I phases 
are very close, A ^ 1.6 and J ^ 0.1, we can drive a CE- 
CO-I to FM-M transition by applying a small magnetic 
field. In the present work we restrict ourselves to two 
choices, J = 0.10 and J = 0.12. While J = 0.1 allows 
closer agreement of temperature and magnetic scales to 
experiments, as seen in Fig. 2. (a), it does not allow much 
room to explore the A dependence. The available A win- 
dow (AA) is ~ 0.1. At the lower end one hits the AF-M 
phase and at the upper end (above the dashed red Hue) 
the CO-state cannot be destabilised using a magnetic 
field. For exploring the BW dependence in more detail 
we choose J = 0.12 allowing a window AA « 0.2. 

We now describe the classification of the regions within 
the zero field CE-CO-I parameter space in terms on their 
responses to field sweeps and rationalize them in terms 
of the zero field phases. This will also allow easier organ- 
isation of the later data. 

The dashed/dotted lines within the CE-CO-I phase de- 
marcate parameter regimes with qualitatively different 
response to magnetic field sweep. The distinct regimes 
are determined by tracking the q — (0, 0) component of 
the magnetic structure factor (5(q)), the Vco and the 
q = (tt, tt) component of _DQ(q) as a function of h. 

(i) The big-dashed line shows the boundary above 
which the CO state does not melt however large the mag- 
netic field. All melting phenomena are confined to the 
CE-CO-I regions below this boundary. Below this line 
there are CE-CO-I regions that either share first order 
boundary with FM-CO for A - 1.6 - 1.65 or with FM- 
M/A-2D at lower A. Melting of the CO can happen by 
a transition to the FM-M phase on applying a field for 
regions close to the FM-M. It is not expected that the 
region A ~ 1.6 — 1.65 which shares boundary with FM- 
CO, will lose CO with applied field. However, we shall 
later see that at intermediate fields the system does phase 
separate between off half doping phases and thereby can 
lose CO. However in the limiting case of very large field 
(where all the spins are polarized), J becomes irrelevant, 
and the half doping FM-CO is recovered back. 



(ii) The finely-dashed line separates regions with dif- 
ferent kinds of hysteretic response. Between the fine and 
the big-dashed lines the CO state melts with increasing 
h and recovers when h is swept back. Between the finely- 
dashed line and the AF(A-2D) boundary the melted CO 
state does not recover when the magnetic field is swept 
back to zero. This non-recovery suggests the existence 
of a metastable FM-M minimum, close by in energy to 
the CE-CO-I absolute minimum even when h — 0. Thus, 
when the field is swept back to zero the system can have 
FM-M volume fractions in the final state. Deeper into 
the CE-region, beyond the finely-dashed line, the FM-M 
minima is no longer metastable at /i = and the CE- 
CO-I is recovered on field reduction. 

(iii) Finally, the dotted line is the boundary between 
inhomogeneous (or partial) melting and homogeneous 
melting. The system melts homogeneously below this 
line while undergoes inhomogeneous melting above it. 
The inhomogeneous melting at fields beyond forward 
melting field, {h^), cannot be rationalized in terms of 
the zero field half doping phases and as mentioned in (i) 
above, is related to phase separation tendency into off 
half doping phases at intermediate magnetic fields. This 
tendency is not just confined to regions with 'FM-CO' 
'CE-CO-I' boundary as mentioned (i), but extends to 
lower Aa as well, where the boundary is between FM-M 
and CE-CO-I . 



B. A - T phase diagrams at J = 0.10 and 0.12 

While the X — J phase diagram gives valuable infor- 
mation on the low T phases, issues like melting temper- 
atures and the thermal evolution of these low T phases 
requires constructing A — T phase diagrams. Together 
with the A — J phase diagram they complete the study 
of the 'clean' system at zero field. 

Here we present two such h = phase diagrams for 
J = 0.10 and J = 0.12. These will provide the Tco 
scales, an estimate of the stability of the CO phase with 
A and will later help to demonstrate sensitivity of the 
field melting scales on small changes in J. Fig. 2(b) and 
(c) show the clean A — T phase diagram, at J = 0.1 and 
J = 0.12 respectively. As a generic feature in both cases, 
for low A, the system has a FM-M ground state, which 
gives way to a relatively narrow A-type AF phase and fi- 
nally to a CE-CO-I state at higher A. Owing to larger AF 
coupling TcE for the J = 0.12 case is larger than for the 
J = 0.10, while the Tc scales are generally suppressed as 
compared to the J = 0.10 case. Starting from the CE- 
CO-I phase at low T, the system goes through a PM-CO 
to a PM-CD phase, with increase in temperature. The in- 
crease in Tco with A signals that beyond a certain A, the 
CO state can be stabilized without CE magnetic order 
and that the stability of the CO state grows with increas- 
ing A. However, at very large A, the electron-phonon cou- 
pling would strongly localise electrons (onsite) and inter- 
site correlations would be weak, suppressing Tco scales. 
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FIG. 3: Colour online: Distinct response to field cycling as a function of A. Top panel shows Vco and ^^(O, 0) and bottom panel 
shows the corresponding resistivity as a function of the applied magnetic field. The X/t values for the systems are indicated and 
Jaf = 0.12. All results are obtained at T=0.02. In the top panel, with increasing A, the forward switching fields /ij increase 
and for A — 1.7, the CO doenot melt, the corresponding resistivities, shown a concurrent switching with the abrupt decrease 
in the Vco and vice verse. Hysteresis is seen in all four cases, except the hysteresis for A = 1.7 occurs only in the magnetic 
sector, while it occurs for Vco, Sfm and p, for the first three CEises. Also note that with increasing A, the melted state has 
successively large amounts of residual CO, which increases to unity for A = 1.7. Further, at intermediate A ~ 1.45, the system 
doenot recover the CO state when the magnetic field is cycled back to zero, signalling that the FM-M is is metastable even at 
zero field. Note, for the A = 1.4 case, because of the structure of the AF(A-2D) phase in Fig. 2(a) one has to go slightly deep 
into the CE-CO phase, where at zero field the CO is recovered in field cycle and hence the FM-M is unstable aX h = 0. 



We do not explore that large A regime. We shall later 
use the trends and the numbers in these phase diagrams 
to make detailed comparison with experiments. 

Having discussed the low T, zero field phase diagrams, 
their thermal evolution and organisation of the parame- 
ter space in terms of the different response to field sweeps, 
we turn to the main theme of the paper which is to de- 
scribe extensively the phenomenon field induced melting 
of the CO state. 



VI. RESULTS AT FINITE FIELD 

For studying the magnetic field effects we track the 
various indicators described in Sec. IV. D, by first cooling 
the system at ft, = and then cycling the magnetic field. 

This section is organized as follows. Sub-sections A 
and B aim at a systematic study of the evolution of the 
field response with A. While in Sub. A we contrast the 
melting vs non melting responses, and through indica- 
tors such as DOS and P(Q) (that measures the Q-Q cor- 
relations on the lattice) hint towards an inhomogeneous 
melted state, in Sub. B, we discuss the complete evo- 
lution with A. In Sub. C we present our h — T phase 
diagrams, explicitly show real space data on inhomoge- 
neous melting in Sub. D and in Sub. E compare our 
results with experiments. 



A. Result of a typical field sweep 

This and the next subsections revolve around Fig. 3. 
This illustrates the field response for four A, J combina- 
tions, with the CO growing stronger as we move from 
(a)-(d). We show the CO volume fraction, Vco and the 
q = (0, 0) spin structure factor in the top panels, and the 
corresponding resistivity p(T, h) in the bottom panels. 

In this subsection we will focus on contrasting between 
the melting vs non melting response. Let us start by il- 
lustrating two CO states (i) at A = 1.55 where the CO 
'melts' is response to increasing h, and (ii) A = 1.70, 
where it does not. These are the last two columns in 
Fig. 3. These two parameter points lie on opposite sides 
of A ~ 1.65 in the A— J phase diagram (Fig. 2(a)). Above 
A ~ 1.65, the FM-CO state does not melt with any ap- 
plied magnetic field, i. e the charge order is not dependent 
on CE magnetic order for its survival. For lower A values 
the CO melts. In both these cases the applied field trans- 
forms the CE order to FM. For the A = 1.55, there is a 
sharp reduction in the CO volume fraction Vco, while 
for A = 1.70 the Vco does not change (there is, however, 
a decrease in the charge disproportionation) . Fig. 3(g)- 
(h) show p{h) at low T. For A = 1.55 there is a large 
drop in resistivity associated with the magnetic transi- 
tion (shown in dashed red lines). The Vco, although 
drops sharply, nevertheless has a finite value. This state 
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FIG. 4: Colour online: Comparison of field response of two 
systems on either side of the critical A('~ 1.65), beyond which 
the CO state doesnot melt with magnetic field, (a)-(b), P{Q) 
vs Q at h/t — and h/t = 0.1 respectively, (c)-(d), N{uj) 
YS ui — fi a.t h/t = and h/t = 0.1 respectively. In (b), the 
dashed line is the P{Q) for A = 1.0, where the ground state 
is a uniform metal. 



with a finite Vco and a 'low' resistivity is likely to be 
a percolative metal. However, this by it self does not 
guarantee that this is the correct ground state, as low T 
field sweeps with multiple competing states can easily get 
trapped into non equilibrium states. While here we take 
this as an indication of possible inhomogeneous melting, 
we carefully separate out cases of equilibrium coexistence 
from nonequilibrium coexistence in later sections before 
making any such claims. Contrast this with the system 
at A = 1.70 where the resistivity continues to be very 
large and the Vco stays unity even at large h. 

For these two cases we now look at distribution P{Q) 
of lattice distortions, and the DOS {N{llj)) for possible 
insight into the melted state. Fig.4.(a)-(b) show P{Q) at 

= and 0.10 respectively. For A = 1.55, at /i = 0, P{Q) 
shows a bimodal distribution. The two values (peaks 
in P{Q)) of distortion on the lattice signifies a possible 
checkerboard CO ground state. This state melts into a 
metal at /i - 0.07, see Fig.3.(g). The P{Q) at h = 0.1, 
Fig. 4(b), shows a broad hump signalling a finite regions 
with localized charges. For comparison, in the same panel 
we have shown data corresponding to A = 1.0 (the dashed 
line), where the ground state is a uniform metal and the 
P{Q) is distinctly peaked close to zero. For A = 1.7, the 
data in panel Fig.4(a)-(b) shows that P{Q) is virtually 
independent of the magnetic field: the high field state is 
a ferromagnetic charge ordered polaronic insulator. 

Fig.4(c)-(d) show the DOS for these A and h combina- 



tions. These correlate well with the P{Q), with the field 
generating a finite DOS at the Fermi level for A = 1.55, 
while there is only a gap reduction (but no closure) for 
A = 1.7. The gap reduction is related to the increase 
in effective BW in the CO phase as the magnetic order 
changes from CE to FM, removing the magnetic blocking 
of electron hopping. 

To sum up: (i) we clearly see two different kinds of 
response to the applied field for systems with different A 
(or BW), and (ii) the Vco and P{Q) results are indicative 
of inhomogeneous melting at intermediate A values. In 
the next subsection we address the A dependence more 
systematically. 



B. The distinct response regimes 

We now take up Fig.3(a)-(b) and their resistivities 
shown in Fig 3(e)-(f), and discuss how they connect with 
the remaining plots in Fig. 3 that we discussed in the last 
sub-section. This will provide the generic nature of re- 
sponse to field sweeping with changing A or BW. 

Fig 3(a) & (e) are at J ~ 0.12 to avoid hitting the 
FM-M phase (see Fig. 2(a)). All the rest are at J = 0.12. 
(i) The data at low A(= 1-40), is shown to illustrate ho- 
mogeneous melting of the CO state as opposed to inho- 
mogeneous melting at higher A. As Fig. 3. (a) shows, the 
Vco reduces to less than 10% for h/t > 0.09. but the 
CO state recovers on reducing h. Note that recovery of 
the CO happens for this 'low' A system even though it 
does not for A ~ 1.45 Fig. 3(b). This is so because 
as seen from Fig. 2(a), at A ~ 1.4 the A-2D phase is 
much wider essentially rendering the FM-M unstable at 
h ^ 0. (ii) For A = 1.45, the A-2D region is small and the 
FM-M remains metastable at /i = and the CO state is 
not recovered even when the field is swept back to zero. 
Further also in Fig.3.(b), the system retains ~ 25% of 
CO volume beyond /ij. (iii) In Fig.3.(c) the system re- 
tains ^ 30% of CO volume beyond and recovers the 
CE-CO-1 state on field reduction, (iv) At even higher 
A ~ 1.7, as discussed earlier, the system retains a global 
CO at all h. Thus, switching field he essentially diverges. 

Let us summarize our findings here. The increase in 
the melting fields with A and systematics of recovery/non 
recovery of the CO with A bear out what was discussed in 
Sec V.A, i.e., if the system is (close to) far from the phase 
boundary (Fig.2(a)), it would (not) be able to recover the 
CO state. And finally with increase in A from 1.4 to 1.65, 
there is an observed monotonic growth in the residual 
CO volume fraction in the field melted state, signalling 
towards inhomogeneous melting. Also if one takes into 
account the resistivity shown in the lower panel, for a 
certain A regime, it would imply that the system is a 
percolative metal. 

While this is indicative of the fact that by appropri- 
ately varying BW one can tune between inhomogeneous 
and homogeneous melting, we would caution the reader 
that this does not amount to saying that the melting 
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FIG. 5: Colour online: The h — T phase diagrams obtained at various A as indicated. At any temperature, the shaded 
regions imply hysteresis, light checkerboard regions imply phase separation and colored areas indicate equilibrium phases. The 
composition of the PS states are demarcated in italics. At low T, the gradual shift of the hysteresis window towards higher fields 
is expected as the CO correlations grow stronger with A. At low T, both at low A(= 1.4) and high A(= 1.7), the intermediate 
field state (beyond the magnetic transition) are uniform states, while for A — 1.4 it is a FM-M, for A — 1.7 it is FM-CO. 
The intermediate A(~ 1.5 — 1.6), the system phase separates beyond ftj, into constituents indicated (in italics). Further, the 
hysteresis window for A = 1.5 extends all the way up to zero fields, implying non recovery of the CO state, increasing A to 
~ 1.6, recovers the CO state. The detailed finite T evolution is considered in the text. 



is inhomogeneous for this entire A(1.4 — 1.65) window. 
Wide domain of metastability of competing phases can 
itself lead to trapping and compHcate genuine phase sepa- 
ration tendencies. This numerical issue needs addressing 
and one needs to carefully disentangle the real inhomoge- 
neous melting (composed of equilibrium phase separated 
states at intermediate fields) from spurious ones, where 
the system can get stuck partly in different metastable 
states leading to nonequilibrium coexistence. We shall 
elaborate upon this in section VII, where we perform al- 
ternate calculations using different protocols to settle this 
issue. 



C. The h — T phase diagrams 

This subsection complements the previous two, in that 
here we provide the thermal evolution of the low T phases 
by constructing the h — T phase diagram. This also al- 
lows for a comparison of trends between theory and the 
available experimental data in the low aA (Ca based) 
manganites as shown in Fig. 1. We will make qualitative 
comparisons in this subsection that will be supported by 
spatial map for the real space evolution of the melting 
process, in Sub. D, and make quantitative comparisons 
are made in Sub. E. 

Fig. 5 shows the h ~ T phase diagrams obtained at 
J = 0.12 and A values indicated. These phase diagrams 
are obtained by sweeping in magnetic fields at various 
temperatures, by first cooling to the relevant temperature 
in zero field. Note, here the definitions of /i+ {h~), used 
to denote the boundary of hysteresis, is taken to be the 
field at which there is a sharp increase (drop) in the FM 
structure factor. This is so chosen because, while for 



Fig.5(a)-(c) these is a CO to CD transition, in 5.(d) there 
is no CO melting and the magnetic transition is a good 
indicator of the switching fields. 

The low temperature phases: The evolution of the 
low T phases with increasing A are exactly the same 
as discussed in Sub. A and Sub.B. The increase in the 
thermodynamic melting fields with increasing A, the 
homogeneous melting and subsequent recovery at low 
A ~ 1.4, the inhomogeneous melting at larger A values 
along with the (non recovery) recovery at intermediate 
(A 1.45-1.54) A ~ 1.55-1.65 and the refusal to melt- 
ing for A > 1.65 are all seen in the h~T phase diagrams 
at the representative A values in Fig. 5. 

Thermal evolution at low fields: An increase in tem- 
perature allows the system to access larger fluctuations 
and thus the hysteresis window narrows with increasing 
temperature. This window tapers off and closes at Tcoi 
with the Tco 's themselves growing progressively with in- 
creasing A, as is also seen from Fig. 2(c). The increase 
in Tco is due to progressive increase in the stability of 
the CO state with increasing A. This increase in stability 
of the CO state also manifests itself in the intermedi- 
ate temperature physics, where at low A ~ 1.4 the loss 
of the CE pattern drives the system metallic, with no 
residual CO correlations. With increasing A however, the 
systems have CO correlations surviving at progressively 
higher temperatures although the long range (q = (tt, tt)) 
correlations get suppressed. The accompanying magnetic 
state at intermediate temperatures is a line like AF phase 
with q = (0,7r) and q = (tt, 0) correlations. This can 
be looked at as a precursor to the low temperature CE 
phase, where the q = (7r/2,7r/2) correlations (crucial 
to the CE phase) do not form at intermediate temper- 
atures. Further increasing the temperature essentially 



10 




FIG. 6: Colour online: The spatial snapshots of charge ordered region (top panel) at A = 1.55, J — 0.12 and the corresponding 
magnetic bonds (lower panel) at various magnetic fields. The red-green checker-board regions are the CO regions while the 
grey regions are the metallic (M) regions. For the lower panel, the red are the FM bonds while the green signifies AF bonds. 
The magnetic field sequence from left to right are h/t — 0, 0.05, 0.06, 0.08, 0.20, 0.08, 0.02, where, first the field is swept 
up starting at zero field, up to a maximum of 0.2 and then swept down back to zero field. Clearly, system recovers the CE-CO-I 
state after going through a percolative metallic phase at intermediate magnetic fields. 



makes the system lose all spin-spin correlations and the 
systems end up in a PM phase. This is accompanied, at 
smaller A 1.4 — 1.5 values, by a metallic state, and at 
higher A ^ 1.6 — 1.7, by a evolution into a metal only 
after going through a CO state. The thermal PM-CO 
to PM-M transition, for this higher A range, signals ther- 
mal liberation of carriers trapped in polaronic state. This 
clearly brings out the stability of CO order at larger A 
quite independent of the spin order. 

Thermal evolution at large fields: Let us contrast the 
low field evolution with that at large finite fields. We 
will later join the two regimes. The expectation that, at 
fields beyond the system will lose all A-2D/CE cor- 
relations, is borne out well for small and large A cases. 
As discussed, while in the small A('^ 1.4) case the CO 
is lost completely at high fields, in the large A(~ 1.7) 
case where the CO does not melt, the systems evolve 
into a FM-CO phase. Increasing temperature will con- 
vert the FM-CO to a FM-M with thermal liberation of 
polaronic entrapments. To the temperatures accessed, 
the high magnetic field is strong enough to resist any 
further thermal transition to a PM state. In contrast to 
the simple evolution at very small (a) and very large(d) 
A, the intermediate A situation (b), (c) is complicated 
by phase separation tendencies. These systems at fields 
beyond phase separate into phases that are at elec- 
tron densities off half doping. The constituents vary with 
changing A, in that, at A ~ 1.5, the coexisting phases are 
'FM-M (ni) + AF-CD (na)' and those at A ~ 1.6 are 
'FM-CO (ni)+AF-CD(n2)'. The densities indicated in 
the brackets imply off half doping densities, typical val- 
ues are given in Sec. VII. A. This phase separation, is the 
origin of inhomogeneous melting and renders the melting 
to b e a continuous transition. The way this tuning can 
happen is discussed in Section VII. These regions are 



indicated, in Fig. 5 (b) and (c), by light colored checker- 
board regions and the phases are indicated in italics. For 
both these A's the corresponding phase separated states 
evolve into a n = 0.5 FM-M. (n=0.5). If one were to 
consider the true limit of large enough field where all the 
spins were fully polarized, then for A < 1.6, the system 
would be a n = 0.5 FM-M and above that it would be a 
n — 0.5 FM-CO (which at high temperature will evolve 
into an FM-M). We relate the intermediate PS to the 
uniform asymptotic limit in Sec. VII. 

Intermediate field thermal evolution: For all the cases 
shown here, the intermediate field shows a transition 
from the respective CE/A-2D/PM to FM state, albeit, 
with the hysteresis region in between, for T < Tqo- 
While for the low A ~ 1.4 and large A 1.7 cases, the 
transitions are benign, for intermediate A ~ 1.5 there is a 
transition at large T and small h, between a PM-M and 
the phase separated state and between the PM-CO and 
the corresponding phase separated state for A 1.6. 

To sum up, the h — T phase diagrams are governed 
by the following principles: (i) The hysteretic nature of 
the field induced transitions, (ii) the increased stability 
of CO correlations with increasing A, (iii) The ensuing 
decoupling of the spin and charge sectors leading to in- 
dependent ordering with increasing A and finally (iv) the 
phase separation tendencies into off half doping phases 
at fields above h^. This of course leaves us to clarify the 
fate of the theoretical spin polarized limit (of h oo), 
where the phase is either CO or metallic, as will be done 
later in the paper. 

To make a few qualitative comparisons with experi- 
mental results shown in Fig. 1, we note that, as for the 
Ca-family, decreasing or equivalently increasing A in- 
creases the CO melting fields and temperatures, as Fur- 
ther, the detailed temperature evolution of low T phases 
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FIG. 7: Colour online: The energies of various phases ob- 
tained as a function of A, measured from the energy of the 
FM-M phase closest to the FM-M, A-AF phase boundary of 
Fig. 2 (a). The black dots indicate the the energy difference 
between points closest to the A-AF region from both the FM- 
M and the CE-CO-I sides, (a) is for J = 0.1 and (b) is for 
J = 0.12. Note here A-AF refers to the hnc-like A-2D phase 
separating the FM-M and the CE-CO-I phases. Also 'closest' 
to a phase, say FM-M refers the closest we could approach 
the phase boundary within numerical accuracy. 



in the h — T phase diagrams in Fig. 5, that have not been 
reported experimentally yet, provides a full picture of the 
'clean' field melting problem. 

Before we make quantitative comparisons with exper- 
iments (in subsection E), we present real space snap- 
shots of the field melting phenomena, which among other 
things lend support to inhomogeneous melting. We turn 
to this in the next subsection. 



D. Spatial character of the melting process 

The spatial data that we provide here are the direct 
snapshots of the real space charge density field and the 
magnetic bonds. These are tracked as the system is swept 
in a magnetic field at low temperature. 

We start with a CE-CO-I system at A(= 1.55). Prom 
the discussion in the previous subsections, we know that 
a system at this A value melts and recovers CO on field 
sweeping, however, the melting appears to be inhomo- 
geneous. The validity that is a genuine inhomogeneous 
melting is provided in Sec. VII. For now we shall assume 
the same and proceed. The spatial patterns are shown in 
Fig. 6, at various fields as indicated in the caption. The 
results are for a 40^ system and have been obtained from 
a run in which h is increased from to 0.2 in steps of 0.01, 
at T = 0.02, and then reduced to zero in the same se- 
quence. The top panel shows the spatial charge ordering, 
while the lower panel shows the corresponding magnetic 
bonds (see figure caption for color convention). 

The primary steps in the evolution are (i) nucleation 
of the FM-M within the CE-CO-I at low fields (as seen 
in the second and third columns from left) (ii) the sharp 



fall in the CO volume fraction, where the system becomes 
a patchwork of FM-CO, AF-CD and FM-M (fourth col- 
umn) (iii) the very large field where the system is com- 
posed of FM-M and FM-CO regions (fifth column) and 
(iv) the subsequent recovery (the rest of the columns). 
While, the nucleation, typical to such transitions, is ex- 
pected, it is interesting to note that the CE phase apart 
from giving way to FM patches, also breaks up into re- 
gions of opposite 'handed' stairs. Since, there is no pre- 
ferred handedness, both being degenerate, small destabil- 
isation of the CE phase results into the system creating 
domains of opposite handed stairs. Within each domain 
the CO state survives, however, its stability is locally 
reduced at these magnetic domain boundaries. This re- 
duction is small and is not apparent in columns two and 
three, where the system started out in the CO state at 
h = 0, however is seen in the last two panels, where the 
system attempts to recover global CO from the melted 
state. Clearly some of the domain boundaries between 
in the CO regions coincide with that of the magnetic do- 
mains boundaries of opposite handed CE regions. This 
fact also makes it rather diSicult to recover the CO com- 
pletely on such a large systems. In column four, the 
system is seen to exist in a coexistent state with FM- 
M, AF-CD and FM-CO coexisting simultaneously. This 
certainly shows explicitly the inhomogeneous nature of 
melting arising out of the phase separation. The column 
five is for h = 0.2 which is large enough to polarize all the 
spins. This is deliberately shown to bring out the effects 
of low lying metastable states that can affect the numer- 
ics. For this value of A the spin polarized problem has 
a global FM-M ground state. However as seen, we find 
a coexisting FM-M and FM-CO. This happens because 
of large domain of metastability of the FM-CO phase, in 
which the system partially gets trapped. This effects the 
results even in column four and one needs an alternate 
calculation to extract the actual constituents of the equi- 
librium phase separated state. This is achieved by doing 
II ~ n calculations, which is the subject of the next sec- 
tion. With hindsight, we mention here, that the correct 
constituents of the phase separation at this A are AF-CD 
and FM-M, both of which are individually at electronic 
densities that are off half doping, as was mentioned in 
subsection C. 

Let us conclude this subsection by briefly summariz- 
ing what happens in such low T field sweeps at other 
A values. For a system with A ^ 1.6, the overall evolu- 
tion would essentially be the same, expect that the phase 
separation is between an AF-CD and an FM-CO, both of 
which are again off half doping phases. For even larger A 
values (> 1.65); the system directly goes from a (n = 0.5) 
CE-CO-I state to a (n = 0.5) FM-CO state. Thus gener- 
ally, for a certain A, J window, the systems undergo in- 
homogeneous melting at intermediate applied magnetic 
fields, where they split into regions of phases of different 
electron densities. This phase coexistent state, can be a 
percolative metal if one of the constituents is an FM-M 
and it manages to percolate. 
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E. Field melting scales 
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We end this section by making a quantitative compar- 
ison of various melting temperatures and melting fields 
with those found in the experiments. The above table 
give the numerical and corresponding values in Kelvin 
of Tco, (hco)' for J/t = 0-1 and J/t = 0.12 for simi- 
lar A values. The numerical Tco values are converted to 
Kelvin, using the kinetic energy scale, t ~ 0.2eVG3, and 
crude finite 2D to 3D scaling of 3/2. The numerical mag- 
netic fields are converted to Kelvin using g^iBhc/t ~ h, 
again with t = 0.2eV. We find that for J/t = 0.1, the 
melting field are much smaller than the corresponding 
melting temperatures. Further, from Fig.l, as a typi- 
cal low disorder example, Pro.sCoo.sMnOs, has Tco = 



250i^, while /ij^ 



27K and Hqq ~ 17K, giving the 



mean hco = 22K. Thus the smallness of the melting 
fields is better achieved if one tunes parameters close to 
J/t = 0.1 as compared to J/t = 0.12. The reason for 
this is seen in Fig. 2(a), where, at J/t — 0.1, the CE-CO- 
I, and FM-M phases are close by, while for J/t = 0.12, 
the two phases are separated by a bigger region of A- 
type AF phase. Fig. 7(a) and (b) show the difference of 
the ground state energies (per site) between the CE-CO- 
I and A- type AF as a function of X/t, from the FM-M 
phase, which is taken to be the reference phase. SE is de- 

fined as 



N 



N being the system size. It is clearly 
seen that the case with J = 0.1 has the CE-CO-I ground 
state energy closer to the FM-M ground state energy, 
while the minimum possible difference between the two 
for J/t = 0.12, is ~ 0.03. Crudely, this barrier to be 
overcome by the applied field at low T. So the smallest 
melting field achievable for J/t = 0.12 would be ~ 42r, 
for X/t ^ 1.4. Where as in principle if one tunes to J close 
to 0.1 such that the A-type AF region vanishes one can 
get zero melting fields. However, we do not pursue the 
direction of fine tuning parameters and focus on study- 
ing the melting phenomenon and the corresponding X/t 
systematics . We have chosen to work at J ^ 0.12 which 
has a larger window oi X/t dependence for this work as 
discussed in Section V-A. 



VII. /i n CALCULATIONS, NATURE OF 
TRANSITION & AN ORGANIZING 
FRAMEWORK 

While the above results on inhomogeneous melting are 
promising, as we saw, a system can get stuck, at low 
temperatures, in states with wide domain of metastabil- 
ity, especially if there are multiple competing phases close 
in energy, as in the present case. Thus there is a need 
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FIG. 8: Colour online: (a) and (b) show the n — fi curves 
at A = 1.50 and A — 1.6 respectively, at various fields as 
indicated. In both cases at intermediate fields, the system 
passes through phase separated states comprising of (AF-CD 
and FM-M) for A = 1.5 and of (AF-CD and FM-CO) for A = 
1.6. Further at high fields both the system are characterized 
my uniform FM-M and FM-CO in (a) and (b) respectively. 
All results are obtained staying at T = 0.02 and J = 0.12. 



to clearly separate genuine inhomogeneous melting from 
any spurious one. This issue is addressed in the follow- 
ing subsection by performing n—fi calculations. This will 
be an alternate way to conclude about phase separation 
at intermediate fields. Once having established the gen- 
uineness of the inhomogeneous melting, we will provide a 
Landau energy landscape as an organizing framework for 
the entire field melting systematics with changing BW. 



A. Field induced phase separation 

We compute the 'chemical potential'- 'electron density' 
(/i — n) curves at various fixed magnetic fields and A val- 
ues. The ^ — n curves are obtained from low temperature 
fi scans of the system, at fixed applied magnetic fields, in 
a protocol that does not retain the memory of previous fi 
steps during the fi sweeps. While runs without memory 
avoid hysteresis effects, the grand canonical calculation 
allows the system to choose the 'best' possible n, thereby 
allowing the system to get into the correct phase at any 
given field. Moreover, we ensured that the system has an- 
nealed well enough by checking that our results hold up 
to large number (8000) of Monte Carlo steps at each fi. 
We also cross checked these results with those obtained 
from a different protocol of coming down in temperature 
and fixed magnetic fields, for a finely discritized set of fi 
values. This ensures that the results are well annealed 
and free from low temperature Monte Carlo problems. 

As seen in Fig. 8 (a), for CE-CO-I systems close to 
the FM-M phase, A 1.5, the CO is quickly lost (be- 
yond h ^ 0.02). At a slightly higher field, the system 
prefers an FM-M state with n = 0.57 up to a certain fi 
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FIG. 9: Schematic sweep rate dependence of the switching at intermediate couplings (A), (a) Equilibrium evolution of the 
magnetization, 5,(0,0), with h (red, solid line), hth is a notional value at which a first order transition would have occurred 
in the absence of PS. The evolution of Vco (blue curves) for different sweep rates is shown in (b)-(d). (b) Slow sweep, (c) fast 
sweep (our regime, see text), and (d) ultrafast sweep. The text discusses the sweep rates in detail. In (c) the shaded region 
depicts hysteresis. 



and then directly go to an ('A-type' AF)-CD phase at 
n = 0.44 leaving the n = 0.5 to be an unstable density. 
Thus, the fixed n{= 0.5)) problem is phase separated and 
infact a percolative metal. The constituents of the PS are 
the FM-M and the AF-CD phases. This is true for all 
systems with A ~ 1.45 — 1.6, at intermediate fields. The 
situation is different for larger A ~ 1.6 — 1.65. For a typ- 
ical case, A ~ 1.6 (Fig. 8 (b)), the system prefers either 
a FM-CO at n 0.52 up to a certain fi and then an AF- 
CD at n — 0.44 at intermediate fields. Such a system 
if forced to be at n = 0.5, as was done previously, will 
phase separate into the above constituents again lead- 
ing to an inhomogeneous melted state. At larger fields, 
both in Fig. 8 (a) and (b), the n = 0.5 state becomes sta- 
ble, recovering the correct asymptotic limits of FM-M for 
A — 1.5 and FM-CO for larger A ~ 1.6. Thus apart from 
ratifying the earlier conclusion of inhomogeneous melt- 
ing, this calculation helps identify the real components 
of the phase separation. Thus, if we look at the spatial 
snapshots for h = 0.08 and h = 0.2 in the field increasing 
run in Fig. 6, we find that the existence of the FM-CO in 
the ground state (for a A = 1.55 system) was purely due 
to it having a wide domain of metastability. This fact is 
not a surprise because, at least at large fields, where the 
J becomes negligible, the FM-M and FM-CO share a first 
order boundary, as seen for small J in Fig. 2 (a) around 
A = 1.6. For the same reason FM-M fractions are trapped 
in the large field FM-CO ground state, for A ~ 1.6—1.65. 
In the fixed density problem, for A < 1.6, we have explic- 
itly checked that annealing the PS system (with FM-M 
AF-CD and FM-CO) to higher temperature first makes 
the FM-CO component disappear rapidly, showing that 
it is purely due to trapping in metastable FM-CO state. 
With this background we now discuss the nature of the 
melting phenomena and its relation to the sweep rate in 
the Monte Carlo and in the experiments. 



B. Nature of field induced transition &: sweep rate 
dependence 

From the — n calculations, it is apparent that the 
melting for a window of couplings is inhomogeneous. 
While for low A ~ 1.4 and high A > 1.65, the respec- 
tive transitions are abrupt, for the intermediate A(~ 
1.45 — 1.65), the melting is inhomogeneous and the transi- 
tion is continuous. Here we discuss a schematic of such a 
transition. We use the ferromagnetic structure factor as 
the order parameter for this discussion. For intermediate 
coupling the magnetic phase separation is between an FM 
and AF states. From their densities one can work out the 
volume fractions of the two constituent magnetic phases. 
Fig. 9(a) shows, the magnetization with increasing ap- 
plied field. The dashed line shows the notional abrupt 
(first order) transition. The continuous line depicts the 
actual transition whereby the magnetization grows con- 
tinuously (from CE-type AF state) with increasing mag- 
netic field to a FM state. The hysteresis occurs in the 
background of this unusual equilibrium physics. Since 
the 'switching' in hysteresis depends on the sweep rate 
let us clarify the experimental and simulation timescales. 

The local relaxation time tioc in electronic systems is 
~ lO"'^^ seconds, but collective relaxation times TcoU, 
say, can be macroscopic, ~ 100 seconds in the CO 
manganites'^. This measurement was at ~ O.QTco and 
TcoU is likely to be much greater at low T. The field cy- 
cling periods Tp^r that we could infer from field melting 
experiments were ^ lOm^^fi!^ Overall r/oc ^ Tper ^ TcoII- 
Our MC results are broadly in the same window. The 
'microscopic' timescale is the MC step. The sweep peri- 
ods were 10'^ — 10^ MC steps (bigger in smaller systems) 
but still ^ 10^^ that one would need to avoid trapping 
in a metastable state. 

The sweep rate dependence of the switching is illus- 
trated schematically in Fig. 9, for an intermediate cou- 
pling system. The left panel, (a), is for a quasistatic 
sweep, Tper S> Tcoii- In this case there would be only pro- 
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FIG. 10: Colour online: The low temperature Landau free energy landscape with the magnetic field axis going into the plane of 
the paper and the other axis, being the general order parameter axis. The three panels (a), (b) and (c) are representative of 
A = 1.5, A = 1.6 and A > 1.65 respectively. All the phases are at electron density n — 0.5, except the ones for which a density 
are shown in brackets along with the name of the phase. The intervening A-type AF region is not shown to avoid cluttering. 



gressive melting and no hysteresis, the system is always in 
equilibrium. Panel (b) illustrates the regime Tper ~ TcoIi, 
where the sweep rate is still 'slow' but the system cannot 
quite track the equilibrium state. In this case there could 
be successive switching. This regime is also out of com- 
putational reach for the system sizes we use. Panel (c) is 
for our regime tioc <C Tp^r ^ t'co/z- The system switches 
at h'^ on field increase, but not necessarily to the under- 
lying equilibrium state. The magnetization, Vco, etc, 
are determined by the presence of metastable states. For 
h ^ , where the equilibrium state is a homogeneous 
FM (at this A) the low temperature system can still re- 
main trapped in the metastable state. Finally, (d) is for 
an ultrafast sweep, r^er ~ tioc, where the system is un- 
able to respond at all to the changing field. 

As shown in panel (c), for sweep rates typical in the 
experiments and in our calculation, the high field state is 
infiuenced by the equilibrium PS and nearby metastable 
phases. So, in any low T field sweep for A = 1.55 (say), 
the observed state is a mixture of the equilibrium AF-M 
-I- FM-M combination and a metastable FM-CO. Increas- 
ing h converts the AF-M to FM-M but the metastable 
FM-CO fraction (also seen in Fig. 6, at /i = 0.2) can be 
removed only by thermal annealing. At larger coupling, 
A > 1.6, the high field low T state would have stable 
(equilibrium) FM charge order. 

C. Landau landscape for melting 

Over the earlier sections, we made a number of conclu- 
sions regarding the BW or A dependence of the response 
to magnetic field sweeps. By invoking a Landau like land- 
scape, with a relevant number of competing phases, we 
organize the broad systematics of the field response with 
changing BW within a single framework. 



From the previous sections we know that the melting 
can either be homogeneous or inhomogeneous. For A less 
than 1.6, the system can melt the CO simply by lowering 
the energy of the FM-M minima with increase in the ap- 
plied field. The increase in field can either lead to a sim- 
ple first order transition, as happens for low A 1.4 or as 
happens for intermediate A ~ 1.5, lead to a field induced 
phase separation. In either case the loss of CO volume 
fraction is guaranteed. However for A ~ 1.6 — 1.65, the 
FM-CO is the one closest in energy to the CE-CO-I (and 
also the true ground state in the limit of h ^ oo) and 
without the intermediate field phase separation, would 
have simple gone over to the FM-CO phase, as happens 
for A > 1.65. Thus the phase separation is crucial for the 
destabilisation of the CO for this A window. 

With this general understanding, let us discuss the 
Landau landscape shown in Fig. 10. This has three pan- 
els depicting the underlying free energy landcapes with 
increasing mangetic fields for three successively increas- 
ing values of A. 

Small A response: Panel (a) of Fig. 10 corresponds to 
A ~ 1.5, where the CO state melts beyond a critical field 
but does not recover when the field is swept back. In 
the energy landscape, the h = landscape has CE-CO-I 
as the global minima, however the FM-M is metastable. 
This metastable FM-M would be necessary to explain the 
non recovery of the CE-CO-I state when the field is swept 
back to zero, as some part of the system remains stuck 
in the metastable state. At h — hi and up to h2, the 
FM-M minima lowers as expected with increasing field. 
If the A is very small ~ 1.4, this continues leading to 
a first order transition to a homogeneous FM-M. How- 
ever, if A ~ 1.5, at h2, the system phase separates into 
off half doping phase (AF-CD -I- FM-M) phases, these 
two minima are depicted in panel (a) at /12. On further 
increasing the field the system evolves to the asymptotic 
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large field ground state of a n = 0.5 FM-M. The phase 
that is closest in energy to this is the FM-CO, as is seen 
in Fig. 2 (a) at low J and A ^ 1.5. Given that the system 
tends to get trapped in the FM-CO, we depict this state 
to be metastable at large fields. 

Intermediate X response: Panel (b) shows a similar 
landscape for A ^ 1.6 — 1.65), range. Here there are 
a few important differences in contrast to that in (a): 
(i) From Fig. 2 (a) the FM-CO phase is the closest in 
energy to the CE-CO-I phase, that can be accessed by a 
magnetic field, (ii) Since we know that the CE-CO phase 
is recovered when the field is swept back, in the h — 
landscape the FM-CO has to be unstable, as opposed to 
the FM-M being metastable at /i = in (a), (iii) the 
phase separation at intermediate fields is between FM- 
CO and AF-CD as depicted, which are again not half 
doping phases, (iv) Finally at large field values the FM- 
CO minima is the global minima with the FM-M minima 
being metastable, as is seen from Fig. 2 (a) at low J. 

Large A response: This is shown in panel (c). Just 
like the very small A systems, the very large A systems 
have a rather benign evolution. As in (b), the phase 
closest in energy to the CE-CO state is the FM-CO state 
and since the CE-CO state is recovered back when the 
field is swept back to zero, this FM-CO state should be 
unstable at ft, = 0. With increasing the field, the FM-CO 
would gradually lower and finally replace the CE-CO as 
the global minima. At large fields (not shown) the CE- 
CO would become unstable. Note here the CO doesnot 
melt and in this view, if the intermediate field induced 
phase separation didnot occur for A ^ 1.6 — 1.65, CO 
melting would not have been possible. 

To sum up, we have qualitatively described the evo- 
lution of the free energy landscape with magnetic field 
and BW variations. Before concluding, we briefly dis- 
cuss some general numerical issues and the large field 
asymptote. 



VIII. ASYMPTOTIC LIMIT & NUMERICAL 
CHECKS 

A. The large h limit 

All through, we have maintained that the large field 
state is rather difficult to capture in the field sweep pro- 
tocol. The large field or spin polarized limit would yield 
the critical A below which the system is a FM-M and 
above which it is a FM-CO. Since this limit is spin po- 
larized, we are effectively at the J = axis of Fig. 2(a). 
There are two equivalent ways of determining the critical 
A. One is setting J = and annealing the system at zero 
field from high to low temperature. The other is to freeze 
all the classical spins to a fixed direction and then to an- 
neal the orbital/lattice variables. Fig. 11 (a) presents the 
A — T phase diagram of the later approach. The parame- 
ter space, at low temperature, is divided into low A < 1.6 
metallic phase and large A > 1.6 CO-00 phase. The 




FIG. 11: Colour online: (a) The A — T phase diagram for the 
frozen classical spin case, (b) The disproportionation 5„ as a 
function of A. 

Tco increases with increasing X/t as is expected. Let us 
now look at the electronic charge disproportionation (5„ 
as a function oi X/t. This is calculated by averaging the 
modulus of the difference of the onsite charge densities 
from n — 0.5 over all sites (5„ = X]i=i n I^-^ ~ 
This although is a bulk measurement, judging by the 
trend of monotonic growth in (5„ in Fig. 11 (b), it is clear 
that even for X/t < 1.6, there is some residual local CO 
regions. We have checked this explicitly by looking at 
snapshots. While above A — 1.6, the system has global 
CO, the residual CO in the lower A case are due to trap- 
ping in metastable state. Thus, although we have cooled 
the system and there is no low temperature sweep in- 
volved, still there is some amount of phase mixing due 
wide domain of metastability of the FM-CO. Thus it is 
bound to affect the low T field sweeps in the earlier calcu- 
lations. We briefly discuss some general numerical issues 
below regarding stability of the results on annealing time 
and mention some general numerical ways to ascertain if 
a phase coexistent state is an equilibrium state or a re- 
sult of a system remaining stuck partially in a metastable 
state. 



B. Numerical checks 

(i) Although TCA gains much in accessing large sys- 
tem sizes, even with TCA performing MC runs with 10^ 
sweeps on large systems at every temperature is pretty 
demanding computationally on large systems. So to en- 
sure that we have annealed well enough, we repeated a 
few calculations on smaller systems employing standard 
exact diagonalization-Monte-Carlo. This was particu- 
larly done at parameter points where the system phase 
separated in the TCA calculations, we presented the 
same state (obtained using the same protocol) to very 
large annealing (~ 40000 MC relaxation steps.) on these 
small systems. The phase separation remained stable 
even after such long runs. 
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(ii) One way to ascertain if the coexistent state (as 
is our intermediate field phase separated states) is the 
equilibrium ground state would be by diagonalizing the 
system in the background of the each of the competing 
phases separately and comparing their energy with that 
of the coexistent state. Still another way would be to 
do alternate calculation which would make the coexistent 
state naturally plausible. In our case the fi—n calculation 
showed the phase separation tendency at intermediate 
fields which unambiguously pointed out the coexistent 
state is naturally favoured over uniform phases. 

IX. CONCLUSIONS 

We have reported the first controlled results on the 
field melting of charge order in half doped manganites 



using an unbiased Monte Carlo method. In this paper 
we have dealt with the 'clean' case, we have mapped 
put the h — T phase diagram exhibiting both hystere- 
sis and re-entrant features. We found the melting to be 
inhomogeneous at intermediate magnetic fields and this 
stands as a testable predictions in this regard. We have 
made extensive comparison between the thermal and the 
magnetic CO melting scales of our calculations and ex- 
periments. Finally, from our numerical results we have 
established the way in which the free energy landscape 
would evolve with increasing magnetic fields and chang- 
ing A (or ta)- We take up the disorder effects on the CO 
melting in the companion papei^^^l. 
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